Origin of the Charge-Orbital Stripe Structure in Lai.^Ca^MnOs (x = 1/2,2/3) 
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We propose the origin of the charge-ordered stripe structure with the orbital ordering observed 
experimentally in Lai-^Ca^MnOs (x = 1/2,2/3), in which the long-range Coulomb interaction 
plays an essential role. We study a Hubbard model with doubly-degenerate e g orbitals, and treat 
the on-site Coulomb interaction (U) and the nearest-neighbor one (V) with the Hartree-Fock approx- 
imation. Both the charge and orbital ordering structures observed in experiments are reproduced 
in a wide region of the U-V phase diagram determined by the present study. The stability of the 
orbital ordering is also confirmed by the perturbation theory. 



PACS numbers: 71.20.Be, 71.10.-w, 71.10.Fd 
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In some perovskite-type hole-doped manganese oxides 
Ai-jAjMnOa (R : rare earth elements, A : alkaline 
earth elements), the colossal magnetoresistance effect has 
been the subject of intense studies. Recently, charac- 
teristic charge ordering phenomena in these materials 
have also attracted a growing interest. Especially, in 
Lai_2;Ca x Mn03 for x — 1/2 and 2/3, it has been re- 
ported that the charge-orbital stripe (COS) structure oc- 
curs in some periodicities which correspond to the com- 
mensurate concentrations fl|-§]. 

Concerning the pure LaMnOs system (x — 0), the an- 
tiferromagnetic (AF) insulating phase appears M, where 
the orbital ordering accompanied with the Jahn- Teller 
(JT) distortion coexists because of the JT active ion 
Mn+ 3 HJ1] ■ In the same way, it is suggested that the JT 
effect is also the origin of the stripe structure in the sys- 
tems with finite carrier concentrations such as x = 1/2 
and 2/3 IHlJlofl. According to this scenario, however, a 
considerably strong JT effect is necessary to realize the 
insulating COS structure |ic|| . In this sense, it is insuffi- 
cient to ascribe the origin of the stripe structure only to 
the JT effect. 

In the present letter, we show that the COS structure 
observed in the Lai-^Ca^MnOa can be explained only 
by considering the Coulomb interaction between carriers. 
We study a Hubbard model with doubly degenerated or- 
bitals which correspond to e g orbitals on Mn ions. By 
treating the Coulomb interaction with the Hartree-Fock 
approximation, we investigate the stability of the COS 
structure observed in experiments for the x = 1/2 system 
Pig. |(a)) @|@ and the x = 2/3 system (Fig. 0(b)) 
pLl|| , and discuss the another type of the charge order- 
ing found by electron microscopy study reported in Ref. 
H (Fig. 0(c)). We show that the COS structure appears 
for the realistic strength of on-site and nearest-neighbor 
Coulomb interactions. In particular, we emphasize that 
the nearest-neighbor Coulomb interaction is indispens- 
able for the occurrence of the stripe structure with the 
orbital ordering. 

Several authors have studied the effect of the on-site 
Coulomb interaction on the orbital ordering in man- 



ganese oxides previously Jl 2| — 11 5[] . Especially, in the x = 
system, it was shown that the orbital ordering is stabi- 
lized by the on-site Coulomb interaction E3]. However, 
the effect of the long-range Coulomb interaction have not 
been studied enough. 

The charge ordering in La 1 _. E Ca ;r Mn03 systems is the 
ordering of Mn ions with different valences: Mn +3 and 
Mn +4 . Three electrons in t^g orbitals construct the lo- 
calized 5 = 3/2 spin. The strong Hund coupling works 
between the localized ti g spin and the S = 1/2 spin of 
the itinerant e g electron. The difference of valences be- 
tween Mn +3 and Mn +4 corresponds to whether the Mn 
ion has an e g electron or not. 

In the COS phase observed in both x = 1/2 and 2/3 
systems, the COS structures are formed in all the a-b 
planes, and they are stacked along the c-axis without 
misfitting. Thus, for simplicity, we investigate the charge 
configuration in the two-dimensional system which cor- 
responds to the a-b plane. The charge ordering is also 
observed in layered- type manganese oxides fl6}| . 

The Hamiltonian which we consider is expressed as fol- 
lows; 
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where rii a denotes a number operator c| c ic< on the site i. 
Indices a and (3 correspond to two orbitals of e g : 3x 2 —r 2 
and 3y 2 — r 2 symbolized by X and Y, respectively. In the 
first term of the Hamiltonian, t"j is the hopping integral 
between the orbital a on the site i and /3 on j. We de- 
fine hopping integrals between nearest-neighbor sites as 
follows; 
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where a; and y denote unit vectors of cc and y directions, 
respectively (Fig. |2|). We introduce values of t and t' as 
t = t and t' = t / A @. Then, the band width W of 
the free system (U = V = 0) is equal to 5to- Hereafter 
we take to as the energy unit . 

In the insulating phase, the screening effect of the 
Coulomb interaction is expected to be suppressed. Thus, 
not only the on-site Coulomb interaction U but also the 
nearest- neighbor one V becomes important. We assume 
for simplicity that the nearest-neighbor Coulomb inter- 
action is not depend on orbitals. 

The Hamiltonian ([!]) is based on the double exchange 
model |p^| . For the strong Hund coupling Jh (~ leV), it 
is expected that the spin of the electron in the e g band is 
fixed to be parallel to the local t 2g spin at T -c Jh- Thus 
we neglect the spin degeneracy for both the ferromagnetic 
and the paramagnetic states. In other words, we consider 
the spinless fermion system |l9[ ]. The present model can 
be interpreted as a single-band Hubbard model with spin- 
dependent hopping integrals provided the orbital degree 
of freedom is expressed as the pseudo-spin. Below, we 
show that the above simplified model can reproduce the 
COS structure observed in manganese oxides. 

The carrier concentration n is given by 

1 1 N 

a=X,Y t=l 

where N denotes the number of sites. In this letter, we 
treat two cases with carrier concentrations n — 1/4 and 
1/6 which correspond to systems with x = 1/2 and 2/3 
in Lai-^CazMnOa, respectively. Hereafter, we apply 
the Hartree-Fock approximation to the two terms with 
Coulomb interactions U and V in the Hamiltonian ([!]), 
and determine U-V phase diagrams. 

First we show the result for the n = 1/4 case which cor- 
responds to the La^Ca^MnOa system. In the phase 
diagram Fig. ||, we see the COS phase, whose structure is 
schematically displayed in Fig. [j](a) or the inset of Fig. [|, 
spreads out in the wide region. The COS phase is realized 
generally in many x = 1/2 compounds of Ri/nAi / 2 Mn03 
(see Fig. 0(a)), and it is called the CE type except for 
the spin configuration. 

There is the ferro-orbital (FO) phase around V = 
for U > 6, where the orbital ordering is realized but the 

charge ordering does not occur; n i ~ 1 and ' ~ 
for all sites. We comment that this FO phase, which ap- 
pears only in the region U > W{= 5), may be the artifact 
of the mean-field approximation. In reality, in a single- 
band Hubbard model, ferromagnetism is hardly realized 
beyond the mean-field approximation p0| . On the other 
hand, the COS phase is realized for U < W, which means 
the validity of the COS phase beyond the mean-field ap- 
proximation. 



For smaller U, the para-orbital (PO) phase is realized, 
where neither charge ordering nor orbital ordering exist 
: nf — nj — n(: const.) for all sites. For the limited re- 
gion of U ~ W and the small V, D and D' phases exist. 
We do not mention their structures here. Note that FO 
and PO phases are metallic. 

Let us turn to the case with the carrier concentration 
n = 1/6. In the recent neutron diffraction experiment for 
La 1 /TCa 2 /3Mn03 ||, the stripe structure corresponding 
Fig. |l](b) is realized undoubtedly. In the present result, 
the same COS structure is reproduced in the fairly wide 
region of the U-V diagram as shown in Fig. [|. We em- 
phasize that the region of the COS phase includes the 
realistic values of J7(~ 5) and V, and it is expected to 
spread further if we include the long-range Coulomb in- 
teraction beyond the nearest-neighbor one. 

For V ~ U, another types of the charge ordering de- 
noted by B and B' occur. However we do not mention 
these structures since they appear only for unrealistically 
larger values of V. Similarly to the n = 1/4 case, there is 
the PO phase with neither charge nor orbital orderings 
for smaller values of U in the n = 1/6 case. The FO 
phase also exists for U > W and V ^ 1. In the n = 1/6 
case, however, the regions for these phases spread out 
more widely than those for the n = 1/4 case. 

The mechanism of the COS structure for n = 1/4 
displayed in Fig. |l|(a) can be understood by the per- 
turbation treatment as follows. We assume t <C U, V 
and set for simplicity t' = in eq. ([|). For large val- 
ues of U and V, it is natural that carriers order alter- 
natively on lattice sites and each orbital on a site is 
singly occupied. In the configuration shown in Fig. 0(a), 
the ground-state energy is given by E g = —2t 2 /{3V) — 
[2/(9UV 2 ) + 29/(45V^ 3 ) + A/{W 2 {U + W)}]t A + 0(t 6 ). 
If the F-orbital is occupied on the O-site (Fig. |l](a)) in- 
stead of X, the ground state energy is raised by AE g = 
{4/(9UV 2 ) + 8/(A5V 3 )}t 4 + 0(t 6 ), since the energy-gain 
through the exchange processes with the y-orbital elec- 
trons at (A,A')-sites are reduced. For the orbital order- 
ing of the n = 1/6 case, a similar argument also holds 
in the order of t 6 provided the charge ordering is con- 
structed. Thus we understand the reason why the long- 
range Coulomb interaction stabilizes the COS structure 
with the orbital ordering. 

Let us discuss another pattern of the COS structure 
shown in Fig. |l](c), which is realized on the surface of the 
sample in Lai/3Ca2/3Cu03 according to Ref. fl. This 
paired COS does not appear in the present phase dia- 
gram (Fig. ^) in which the JT distortion is neglected. 
Concerning the JT distortion, each non-paired Mn 3+ 
stripe causes significant lattice distortion in its immediate 
neighborhood and the overall strain energy will be low- 
ered by forming periodic array of pairs on Mn 3+ stripes 
separated by undistorted regions of Mn 4+ ions ||. We 
note that the JT distortion is stronger than the bulk one, 
since the elastic constant on the surface is smaller. 

In order to estimate the strength of the JT distortion 
energy enough to construct the paired stripe structure, 
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we consider the following JT Hamiltonian Hjt- 



Wjt = Q 
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n iY ) - (rijx - n jY )}, (6) 



where ((••■}} denotes the sites which correspond to the 
paired stripe shown in Fig. |l|(c); in a pair ({i, j)), sites i 
and j are the left and right site of the pair, respectively. 
In the Hamiltonian (||), Q denotes the potential energy 
which corresponds to the electron-lattice coupling related 
to the linear displacement of the JT distortion ]21J . For 
simplicity, we neglect the mixing term between two or- 



and 3y 



to realize the configuration 



bitals 3x 2 

in Fig. [l](c). The lattice elastic energy, which is propor- 
tional to the square of the displacement, is neglected. 
For the system added the above Hamiltonian (|J), we cal- 
culate the energy under the constraint; riix = n^y and 
niY = rijx for Figure || shows the contour di- 

agram for values of Q which are needed to realize the 
paired stripe structure for certain values of (U,V) with 
the same ground-state energy of the Q = system with- 
out the constraint. 

Obviously, in the case of V <C U, fairly large values of 
Q are needed to realize the paired stripe structure. On 
the other hand, for the region of the COS phase in Fig. 
|], the paired stripe structure can be realized with rela- 
tively small values of Q. This indicates that the nearest- 
neighbor Coulomb interaction is also important to con- 
struct the paired stripe structure. The result suggests 
the possibility that the paired stripe structure reported 
in Ref. j5j on the surface is realized owing to the smaller 
elastic constant on the surface than that in the bulk. 

Finally, we explain some experiments on the COS 
phase and discuss the validity of the long-range Coulomb 
interaction scenario; in Pri-^Ca^MnOa systems, it was 
reported that the charge-ordering transition temperature 
(^co) is suppressed with external pressure |22|j . Un- 
der pressure, the transfer integral is expected to become 
large, and both U/to and V/to decrease. In fact, ab- 
solute values of the resistivity are suppressed with in- 
creasing pressure 22 . The reason why the charge order- 
ing becomes unstable is expected to be the reduction of 
U/to and V/to in the phase diagram, Figs. || and [|. The 
present conclusion is qualitatively consistent with the 
above experimental result. The similar behavior is also 
reported in layered-type manganese oxides with chemical 
pressure [ p3[ . 

On the other hand, recently, the charge-ordering phe- 
nomenon was observed in thin films (thicknesses are 
500 - 2000 A) in Nd .5Sro.5Mn0 3 ||. It was shown 
that Tqo almost docs not change when the thickness of 
the sample changes, although the JT distortion varies 
sensitively with the thickness of the sample. This is con- 
sistent with the present result because we show that the 
COS structure can appear without taking account of the 
JT effect. It is suggested that the JT effect is not so 
important to the construction of the COS structure ob- 
served in manganese oxides with x > 0.5. 



We would like to comment on that we treated the spin- 
less system in the present study. In manganese oxides in 
which the charge-ordered phase is observed, Tqo is higher 
than the AF transition temperature for x > 0.5 ||, 
which will mean that the AF ordering is a phenomenon 
in the lower energy scale. It is of course expected that 
the RKKY interaction, which is neglected in the present 
study, becomes effective to the AF ordering at lower tem- 
peratures. 

In summary, we have shown that the COS phase ap- 
pears in a wide region of the U-V phase diagram for 
both n = 1/4 and 1/6 cases. We conclude that the COS 
structure observed in La^Cai-^MnOa for x = 1/2 and 
2/3 is originated from the Coulomb interaction. It is 
emphasized that doubly-degenerated e g orbitals and the 
nearest-neighbor Coulomb interaction are important to 
realize the COS structure. According to the present sce- 
nario, the COS structure is realized by the long-range 
Coulomb interaction at first. After that, the JT distor- 
tion occurs so that the COS structure is stabilized. In 
conclusion, the JT distortion is the consequence of the 
COS structure but not the origin of it. 
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FIG. 3. U-V phase diagram for the n = 1/4 case. The 
COS phase corresponds to the COS phase with the orbital or- 
dering. Only the FO and PO phases are metallic. The inset 
shows the schematic charge configuration in the COS phase 
(see also Fig. 1(a)). 
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FIG. 4. U-V phase diagram for the n = 1/6 case. Only 
the FO and PO phases are metallic. The inset shows the 
schematic charge configuration in the COS phase for the 
n = 1/6 case (see also Fig. 1(b)). 



FIG. 1. (a) Schematic configuration of the COS structure 
observed in experiments for (a) the x = 1/2 system and (b) 
the x — 2/3 system, (c) Paired COS structure reported in 
Ref. [5] for the x = 2/3 system. 
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FIG. 2. Hopping integrals between nearest-neighboring 
sites. 



FIG. 5. Contour diagram for values of Q which are needed 
to realize the paired COS structure for (U, V). For detail, see 
the text. 
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